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ABSTRACT 


Observations and theories of particle acceleration in solar flares are 
reviewed. The most direct signatures of particle acceleration in flares are 
gamma rays, X rays and radio emissions produced by the energetic particles in 
the solar atmosphere and energetic particles detected in interplanetary space 
and in the Earth's atmosphere. We discuss the implication of these 
observations, we review the theories of stochastic and shock acceleration as 
well as acceleration in direct electric fields, and we briefly discuss 
interplanetary particle propagation. We attempt to present an overview of the 
highlights of both current and promising future research. 
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I. INTRODUCTION 
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Acceleration of energetic particles is a wide spread phenomenon in 
nature, one that occurs at a variety of sites ranging from the Earth's 
magnetosphere to distant objects such as supernovae, active galaxies and 
quasars. There are in fact many explosive phenomena in astrophysics, solar- 
flares among them, in which energetic particles are routinely produced and 
often contain a large fraction of all the available energy. 

It is widely believed (e.g. Syrovatskii 1981) that solar flares draw 
their energy from the annihilation of magnetic fields. The strong electric 
fields which accompany such annihilation should accelerate particles, and in 
addition the rapid energy deposition following the annihilation should be an 
important source of shocks and turbulence. As proposed by Fermi (1949), 
charged particles can be accelerated to high energies in repeated reflections 
from magnetized clouds, or, in the more recent view, from hydromagneti c 
turbulence and shocks. This mechanism must be important in solar flares where 
shocks are known to exist and turbulence is expected to be produced by both 
shocks and other mechanisms. 

Because the Sun is close to the Earth, particle acceleration in flares 
can be observed in considerable detail. relevant observables are 


electromagnetic radiations produced by the ccelerated particles in the solar- 
atmosphere and energetic particles that escape from the Sun and are directly 
detected in interplanetary space. From the analysis of electromagnetic 
radiations and energetic particles observed from flares it follows that not 
all solar flare particles species are accelerated to their final energies at 
the same time. There appears to be a first-phase particle acceleration which 
produces mostly non-relativistic electrons and a second phase which 
accelerates ions and relativistic electrons. It is not clear, however. 
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whether these two phases are manifestations of different acceleration 
mechanisms, or whether they are due to the evolution in space and time of the 
same mechanism modified by different enerqy loss and particle transport 
processes. Also, recent gamma-ray observations have shown that energetic ions 
can be produced very promptly, within a few seconds of the nonrelativistic 
electrons. 

The energy spectrum resulting from solar flare particle acceleration can 
be determined down to a few MeV/nucleon from the maximum flux at each energy 
as observed in interplanetary space. This spectrum is clearly not a power 
law; simple shock and stochastic acceleration models appear to produce eneray 
spectra that can fit the data. The chemical composition of solar energetic 
ions observed in interplanetary space, while grossly similar to the 
composition of the solar atmosphere, varies from one flare to another, and 
occasionally departs dramatically from the photospheric composition. These 
variations probably result from the existence of acceleration thresholds and 
injection mechanisms which depend on the charge and mass of the particles. 

The ionic states of energetic particles give some indication that the 
particles observed in interplanetary space are accelerated in the corona; they 
also seem to imply that selective heating mechanisms are occasionally 
important. 

We review the observations and their implications in Section II, we treat 
acceleration by turbulence, by shocks, and in direct electric fields in 
Section III, we discuss the problem of the determination of solar flare 
particle energy spectra from interplanetary observations in Section IV, and we 
summarize the paper and provide an outlook for future research in Section V. 
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II. ENERGETIC PARTICLES IN SOLAR FLARES 

The solar flare 1s capable of producing a rich and complex particle 
population characterized by a broad spectrum of particle energies and 
containing many different cherical. Isotopic and Ionic species. Of these, the 
most commonly accelerated particles are electrons of energies In the range 
from about 10 to 100 keV. Information on these particles Is obtained from 

I 

hard X-ray observations (e.g. Hudson 1979, Kane et al . 1980) . from microwave | 

and Type III radio observations (e.g. Lin 1974), and from direct electron | 

detections In Interplanetary space (e.g. Lin 1974, Ramaty et al. 1980). j 

i 

Solar flares also accelerate protons, nuclei and relativistic I 

electrons. The nucleonic component Is observed In interplanetary space (e.g. 

McDonald, Fichtel and Fisk 1974), It Is occasionally detected on the ground by 
neutron monitors (e.g. Shea and Smart 1973, Duggal 1979) and Its interactions 
at the Sun are manifest in gamma-ray lines (Chupp et al. 1973, Ramaty, 

Kozlovsky and Lingenfelter 1975, Chupp 1982). Relativistic electrons are 
observed directly in interplanetary space (Simnett 1971, Datlowe 1971, 

Evenson, Meyer and Yanagita 1981) and their interactions at the Sun produce 
gamma-ray continuum (Surl et al . 1975) and Type IV radio emission (e.g. Ramaty 
et al. 1980). 

In the present Section we review the most pertinent data on solar 
energetic particles and discuss their implications. Electromagnetic emissions 
are discussed In Ila and the direct particle observations are in Ilb. Since 
much of the information on solar radio emission is treated in a separate 
article In this volume (Smith and Gtoldman 1982), we discuss this topic only 
briefly here. Solar hard X-ray observations have also been reviewed in 
considerable detail recently (Hudson 1979, Kane et al. 1980) as have been the 
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gamma-ray observations (Ramaty 1982, In this volume). 
a. Electromagnetic Radiations 

The electromagnetic emissions that contain the least ambiguous 
Information on solar energetic particles are radio emissions, X-rays and 
gamma-rays. Other emission such as EUV and white light radiations are not 
discussed In the present article, but the reader can refer to other recent 
reviews (e.g. Kane et al . 1980, Ramaty et al. 1980). 

It has been known for sometime (Wild, Smerd and Weiss 1963, Frost and 
Dennis 1971, Lin lO'^^, Ramaty et al . 1980) that at least two phases of 
particle acceleration can occur In solar flares. The first phase Is 
characterized by Impulsive bursts of hard X-rays and microwaves and by Type 
lii radio emission, while the second phase Is manifest In Type II bursts, 
microwave and metric Type IV emission, flare continuum radio emission and 
gradual hard X-rays. Gamma-ray emission is probably due to particles 
accelerated In both the first and second phases. The first phase has typical 
rise times of seconds or less and durations as short as a few seconds, while 
the second phase. In general, has longer rise times and durations. We proceed 
now to discuss these emissions In more detail. 

Radio Emissions 

As just mentioned, first-phase acceleration In the radio band Is 
characterized by Impulsive microwave bursts and Type III radio bursts. 
Impulsive microwave bursts are generally believed to be due to gyrosynchrotron 
radiation of electrons of energies In the range from several tens to several 
hundreds of keV (Takakura 1960, Ramaty 1969). Type III bursts result from the 
conversion Into electromagnetic emission of Langmuir waves excited by beams of 
electrons of energies In the 10 to 100 keV range (Wild, Smerd amd Weiss 



OrtiGiNAL PAGE IS 
OR POOR QUALITY 


8 


1963). The same energy electrons also produce hard X-rays (see below). These 
electrons have access to both closed field lines of low loops, where they 
radiate microwaves and hard X-rays, and open field lines where they produce 
type III bursts. The electrons on open field lines can also be directly 
observed in interplanetary space (e.g. Lin 1974). 

Second phase acceleration is characterized by Type II radio bursts, 
microwave Type IV bursts, flare continuum and moving Type IV emission. Type 
II bursts, presumably also due to the conversion of Langmuir waves, are 
indicators of shock fronts (Wild, Smerd and Weiss 1963) in the corona. 

Particle acceleration by these shocks, or by the turbulence behind them, is 
consistent with many second-phase phenomena including the protons and nuclei 
observed in interplanetary space. This conclusion is based on the good 
correlation between proton events and Type II bursts (Svestka and Fritzova 
1974). The Type II emission itself is probably produced by 10 to 100 keV 
electrons that are also accelerated by the shock. 

Type IV emission, in general, is due to gyrosynchrotron radiation of 
relativistic electrons. This emission includes microwave Type IV bursts, 
produced in low coronal loops where the magnetic field is relatively high, and 
moving Type IV bursts which are generally observed in the corona. The 
microwave Type IV emission can be distinguished from impulsive microwave 
bursts by their delayed rise to maximum intensity. The moving Type IV bursts 
are due to synchrotron radiation of ~ 1 MeV electrons in coronal magnetic 
fields of a few gauss. Flare continuum emission is also observed in the 
corona at the time of the passage of a Type II burst. As opposed to moving 
Type IV bursts, this emission is stationary and is believed to be due to 
plasma radiation (Dulk, Melrose and Smerd 1978). 
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Hard X-Rays 

Hard X-rays In solar flares result from the bremsstrahlung of electrons 
In the energy range from about 10 to several hundred keV. In the X-ray band, 
first-phase acceleration Is characterized by Impulsive hard X-ray bursts. The 
correlation between these bursts and Impulsive microwave bursts (e.g. Lin 
1974) Indicates that both emissions should be produced by essentially the same 
electron population. As already mentioned, these electrons also produce Type 
III radio bursts. The first-phase mechanism Is thought to be responsible for 
the Impulsive acceleration of 10 to several hundred keV electrons In flares 
(e.g. Ramaty et al. 1980). 

The most detailed Information on first phase acceleration of electrons 
comes from hard X-ray observations. This Information, however. Is model 
dependent. Two limiting cases exist regarding the nature of the electrons. 
They could •'dong lu a n-^nthermal population whose number density Is lower 
than that of the relatively cool ambient medium with which It Interacts, or 
they could form a quasithermal hot plasma. In the former case, the radiation 
yield Is the ratio of the bremsstrahlung production rate to the nonradlative 
colHslonal loss rate. Since In the X-ray band this yield Is very small (' 
10"5 at 25 keV), hard X-ray production In flares should be accompanied by the 
deposition of large amounts of energy Into the solar atmospher-e. On the other- 
hand, In thermal models (Chubb 1970, Crannell et al . 1978), a larger fraction 
of the electron energy content could In principle be emitted as hard X-rays. 
But It turns out that It Is difficult to confine the hot electrons to dense 
regions for times comparable to their radiation loss time (Smith and Auer- 
1980) and therefore the radiation yield of the thermal models is probably not 
much higher- than that of the non-thermal one. 

There are also a variety of particle Interaction models which can be 
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crudely classified as thin- or thick-target {e.g. Karre et al. 1980). In the 
thin-target model (Datlowe and Lin 1973) the X-rays are produced by electrons 
which escape from the Interaction region at the Sun, while In the thick-target 
model (Brown 1971, Hudson 1972), X-ray production takes place as the electrons 
slow down In the solar atmosphere. Clea.ly, the radiation yield Is lower for 
a thin-target than for a thick-target, because for the latter, in addition to 
the colllslonal losses, the electrons also carry away kinetic energy. 

However, from the comparison of number of electrons observed In Interplanetary 
space with that needed to produce the Impulsive hard X-rays at the Sun, It 
follows (Lin 1974) that the majority of the 10 to 100 keV electrons remain 
trapped at the Sun and produce X-rays In thick-target Interactions. 

Various estimates exist of the energies deposited In the solar atmosphere 
by electrons accelerated In the first phase. From the summary of Ramaty et 
al. (1980) we have that In the thick-target and nonthermal model, the energy 
deposited by electrons above 25 keV Is ^ 2xl032 ergs for the August 4, 1972 
flare, which was one of the strongest events observed. The energy deposited 
In other cases ranges from about 5 to 50x1029 ergs In typical strong events, 
and from about 2 to 20xl028 ergs In typical small events. The fact that these 
energies are comparable to the total observed flare energies has led to the 
suggestion (Lfn and Hudson 1976) that most flare phenomena could be accounted 
for, at least energetically , by the Interaction of the first-phase electrons 
with the solar atmosphere. 

First-phase acceleration must be very efficient, both In the amount of 
energy that It converts Into accelerated particle energy, as we have just 
seen, and the number of ambient particles that It accelerates. The mechanism, 
however, accelerates fewer protons to 10 to 100 keV than electrons, as 
Indicated by the absence of nonthermal wings of hydrogen Lya (Canfield and 
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Cook 1978). We discuss below the recent g«nma-r*ay evidence for the Ifnpulslve 
acceleration of MeV protons. As we shall see, the energy content In these 
particles Is also lower than that In the 10 to 100 keV electrons. Mechanisms 
for first-phase acceleration hav.-» been reviewed by Smith (1979). 

In addition to the Impulsive ccxnponent, so'ar hard X-rays occasionally 
exhibit a gradual component (duration ^ 10 minutes) which Is believed to be a 
manifestation of second-phase acceleration. Hudson (1979) has reviewed the 
evidence fur this component. It can be four.1 In hard X-ray observations from 
flares located behind the limb of the Sun (Hudson, Lin and Stewart 1981) and 
from flares on the visible solar hemisphere which produce extended X-ray 
bursts (Frost and Dennis 1971, Hoyng, Brown and Van Seek 1976). These X-rays 
are produced In the corona and their spectrum Is generally harder than that of 
the Impulsive X-rays. This Indicates additional acceleration. But the energy 
deposited in the solar atmosphere by the 10 to 100 keV electrons responsible 
for the gradual X-rays 1$ lower by at least an order of magnitude than that 
deposited by the electrons of the Impulsive phase (Hudson 1979). 

Gamma-Rays 

Gamma-ray lines, produced In energetic particle reactions In the solar- 
atmosphere (Ramaty, Kozlovsky and Lingenfelter 1975), are tracers of the 
nucleonic component accelerated in flares. Many narrow and broad lines are 
produced (e.g. Ramaty, Kozlovsky and Lingenfelter 1979). The strongest narrow 
line Is at 2.22 MeV from -leutron capture on hydrogen. This line has been seen 
from a number of flares so far, as have some of the other strong narrow lines, 
e.g. at 4.44 MeV from at 6.13 MeV from 16o and at 0.51 MeV from positron 
annihilation (Chupp 1982, Ramaty 1982, In this volume). In addition to being 
observed individually, nuclear lines make a significant contribution to the 
total gamma-ray emission above 1 MeV (Ramaty et al. 1980). In particular. In 
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the 4 to 7 MeV band, nuclear radiation appears to be che dominant emission 
mechanism (Ramaty, Kozlovsky and Surl 1977, Ibrf^lmov and Kocharov 1977). On 
the other hand, gamma-ray emission up to an MeV Is almost entirely due to 
electron bremsstrahlung and hence is a good tracer of the relativistic 
electrons In this energy range. Nuclear reactions In the solar atmosphere 
occasionally produce high energy neutrons which can travel as far as the Earth 
resulting In detectable neutron fluxes (Lingenfel ter and Ramaty 1967). High 
energy neutrons were observed from a large flare In 1980 (Chupp et al. 1982). 

Gamma-ray observations can provide Information on the timing of the 
acceleration of protons, nuclei and relativistic electrons, they can determine 
the energy deposited by the nucleonic component In the solar atmosphere, and 
they can place constraints on the energy spectra of the protons and nuclei. 
Through Doppler shifts, gamma-ray lines could provide unique Information on 
the beaming of the energetic particles (Ramaty and Crannell 1976, Kozlovsky 
and Ramaty 1977, Zwelbel and Haher 1982). 

We first address the question of the timing of the nucleonic acceleration 
and Its relationship to the timing of the acc?»lerat1on of th*» 10 to 100 keV 
electrons In the first phase. This question was studied (Bal and Ramaty 1976, 
Lin and Hudson 1976) for the August 4, 1972 flare which was the first flare 
from which nuclear gai.^-rays were seen (Chupp et al . 1973). 

The time dependences o^ hard X-rays, gamma-r jy continuum, high fr'^quency 
microwave emission and the 2.22 MeV line for the August 4, 1972 flare are 
shown In Figure 1 (from Bai and Ramaty 1976). As can be seen, the emissions 
produced by relativistic electrons (> 0.35 MeV gamma-ray continuum and 
high-frequency microwave emission/ had a time history different from that of 
the X-ray emission produced by < 100 keV electrons. This difference Is 
manifest In that the gamma-ray continuum and microwave emission reached peak 
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strength a few minutes later than the X-rays. In addition, as can be seen 
from Figure 1, the time profile of the 2.22 MeV line could be better explained 
when the time profile of neutron production was assumed to be similar to the 
time profile of the gamma-ray continuum rather than that of the X-rays. (See 
article by Ramaty 1982, In this volume, for a detailed discussion on the delay 
of the emission of 2.22 MeV photons owing to the finite capture time of the 
neutrons). This led Bal and Ramaty (1976) to conclude that the >0.35 MeV 
electrons and > 10 MeV protons are accelerated In the second phase. However, 
as can be seen In Figure 1, even though the build-up of the higher energy 
emissions In the August 4, 1972 flare was slower than that of the hard X-rays, 
the two emissions were not clearly separated In time. 

As opposed to the gradual build-up of the nucleonic component on a time 
scale of a few minutes seen In the August 4 event, the recent SMM (Chupp 1982) 
and HEAO-3 (Prince et ol . 1982) observations have provided examples of very 
prompt acceleration of protons and nuclei. The most useful Infor-mation in 
this regard comes from the comparison of the time history of the prompt gamma 
rays In the energy band from 4 to 7 MeV (which Is dominated by nuclear gamma- 
rays) with that of the hard X-rays. For the June 7, 1980 flare (Chupp 1982) 
this comparison (Bal 1982) Implies that the < 100 keV electrons and > 10 MeV 
nuclei were accelerated In very close time proximity (less than a few 
seconds). There are, nevertheless, differences between the two time 
profiles: there Is a delay of approximately 2 sec between Individual hard 

X-ray and gamma-ray peaks, and the gamma rays reach maximum strength at least 
10 seconds after the X-rays. The small but finite lag of gamma rays behind 
the hard X-rays In the June 7 flare would be yet another manifestation of a 
second acceleration phase. This type of second-phase acceleration, I.e. one 
which Is very closely correlated with the first-phase acceleration, has baen 
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referred to as second-step acceleration (Bai and Raraaty 1979, Bal 1982). 

We next discuss the energy content in the nucleonic component and its 
relationship to the energy deposited by electrons accelerated in the first 
phase. The nucleonic energy content is the sum of the energy deposited during 
gamma-ray line production and the energy carried away by particles escaping 
from the Sun. 

We consider numerical values for two flares (August 4, 1972 and June 7, 
1980) for whic.i the gamma-ray data have been analysed in considerable 
detail. These are given in Table 1. Here (>lMeV/nucleon) and Wgsc 
(>lMeV/nucleon) are, respectively, the energies deposited and carried away by 
the nucleonic component, (>20keV) is the energy deposited by electrons 
above 20 keV, and Np (>10 MeV) and (>10 MeV) are the numbers of protons 

above 10 MeV that interact and escape from the Sun, respectively. For the 

August 4 flare, W,j and Np are from Ramaty (1982, in this volume), and 
Ngsc,p are from Lin and Hudson (1976) and We has been given earlier in this 
section, ."or the ^ ,ie 7 lare, Wj and Np are again from Ramaty (1982, in this 

volume), Wgsc and Ngjg^p are from Von Rosenvinge, Ramaty and Reames (1981) and 

Wg is from A. Kiplinger (private communication 1981). 

Considering the numerical values of Table 1, we first note that the 
escape conditions for the nucleonic component can vary considerably from one 
flare to another. While for the June 7 flare only a small fraction of the 
particles escaped from the Sun, for the August 4 flare the number of protons 
observed in space exceeded the number that interacted at the Sun by more than 
an order of magnitude. Evidently, in certain cases particles have access to 
open field lines while in others they do not. In both cases, however, the 
total energy in the nucleonic component is less than the energy deposited by 
the electrons accelerated in the first phase. Shocks, turbulence, or both. 



produced during energy deposition In the first phase, could thus be 
responsible for second-phase acceleration. 

The gamma-ray observations also set constraints on the energy spectra of 
the protons and nuclei. The derived energy spectra (Ramaty 1982, In this 
volume) turn out to be consistent with those obtained from direct particle 
observations In Interplanetary space (McGuire, Von Rosenvinge and McDonald 
1981 and Ilb). This agreement implies that the same acceleration mechanism 
could be responsible for both the escaping and trapped nucleonic component, 
even though the ratio between the numbers of particles In these two 
populations can vary dramatically from flare to flare. On the other hand we 
cannot rule out the possibility that these two populations are produced by 
different mechanisms. 

For a discussion of additional Implications of the solar gamma-ray 
observations (e.g. beaming, photospheric ^He abundance, limb darkening 
effects) the reader is referred to the article by Ramaty (1982, in this 
vol ume. 

b. Energetic Particles 

Complementary to the electromagnetic radiations discussed in the previous 
subsection are the direct particle observations In Interplanetary space and In 
the Earth's atmosphere. Spacecraft observations (e.g. Van Hollebeke 1979, 
Gloeckler 1979) can determine the energy spectra of the various particle 
species Including electrons, the abundances of essentially all elements from H 
through Ml, and Isotopic and lorilc abundances for a few abundant elements. 
Ground based neutron monitor observations (e.g. Duggal 1979) extend the 
spectrum of the nucleonic component of solar flare particles to about 20 
GeV. Particle observations, however, can only provide limited Information on 
the timing of the acceleration. We proceed now to discuss the most pertinent 
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Energy Spectra and Electron-Proton Correlations 

Detailed studies of solar flare proton spectra that attempted to 

eliminate the effects of coronal and interplanetary transport were made by 

Bryant et al. (1965), Reinhard and Wibberenz (1974) and Van Hollebeke, Ma Sung 

and McDonald (1975). In the study of Van Hollebeke et al. the effects of 

propagation on the observed spectra were minimized by considering only 
particle events from flares that were well connected magnetically to the 
observing spacecraft and by constructing the proton energy spectra at times of 
maximum intensity at each energy (see Section IV). Because of instrumental 
limitations, however, the resultant energy spectra extended only over the 
narrow energy range from 20 to 80 MeV and therefore could not differentiate 
between various possible spectral fits such as power laws or exponentials. 

A more recent study (McGuire, Von Rosenvinge and McDonald 1981) has 
extended considerably the particle energy range under investigation (1 to 400 
MeV for protons) and has provided data on a-particle spectra as well. By 
using the same techniques for minimizing the propagation effects as Van 
Hollebeke et al. (1975), McGuire et al. (1981) find two spectral forms that 
provide good fits to the data. These are the Bessel function spectrum (Ramaty 
1979, Section Ilia) 

^(E) cc g K2[2(36/aT)^^] , (1) 


and an exponential spectrum in rigidity 
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Here dJ/cE is differential intensity measured in partial es/(cm2sec sr MeV/ 
nucleon), v = CB, E and R are particle velocity, energy per nucleon and 
rigidity, respectively, and K2 the modified Bessel function of order 2 
(Abramowitz and Stegun 1966). The parameters aT for the Bessel function and 
Rq for the exponential characterize the shape of the particle energy 
spectra. Equation (1) was shown to be the solution of a transport equation 
for stochastic acceleration with acceleration efficiency a and loss time T (M. 
A. Lee private communication 1978, Ramaty 1979, Section Ilia). 

An example of the results of McGuire et al. (1981) is shown in Figure 

2. Here the data points show the observed particle intensities from the 

November 5, 1974 event, and the solid and dashed curves are given by equations 

(1) and (2), respectively. The spectral parameters for this event are aT = 

0.024 and R = 73 MV for protons, and aT = 0.015 and R = 80 MV for a 
0 0 

particles. A very important result of McGuire et al. (1981) is that power 
laws in energy (dJ(E)/dE « e“^ with y * constant) do not fit the data over the 
entire observed energy range, as can be clearly seen from Figure 2. 

The variability of the spectral parameters aT and Rq from one flare to 
another was also studied by McGuire et al. (1981). They find, as did Van 
Hollebeke et al. (1975) before, that for the well-connected flares the 
spectral parameters are confined to rather narrow ranges. For protons, aT = 
0.025 ± 0.011 and Rq = 70 ± 27 MV. The a-parti cl e spectra are generally 
steeper than the proton spectra. The ratios of the proton to a-particle aT's 
and Rq's are 1.6 ± 0.2 and 0.7 ± 0.1, respectively (R. McGuire, private 
communication 1981). 

We pointed out in I la that constraints on the energy spectra of the 
protons and nuclei that interact at the Sun can be set by gamma-ray line ratio 
observations. Unlike the interplanetary observations, these data are not 
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influenced by propagation effects, but they depend on the assumed interaction 
model. For thick-target interactions, Ramaty (1982, in this volume) finds 
that for 7 flares for which line ratios are available, aT (assumed to be the 
same for all particle species) is in the range 0.014 to 0.020. This 
relatively narrow range is consistent with that found by McGuire et al. (1981) 
from interplanetary observations. As already mentioned, the implication of 
this result is that a single mechanism could be responsible for both the 
escaping and trapped nucleonic component. 

We turn now to a discussion of electron energy spectra and correlations 
between electrons and protons of various energies. Energy spectra of 
electrons accelerated in solar flares were analysed by Lin (1971, 1974) and 
recently by Lin, Mewaldt and Van Hollebeke (1982) who have minimized the 
propagation effects in the same fashion as was done for the protons and a 
particles discussed above. An energy spectrum of a large electron event is 
shown in Figure 3 (from Lin, Mewaldt and Van Hollebeke 1982). The break at ^ 
100 keV, characteristic of all solar flare electron spectra, is probably due 
to the acceleration process itself (Lin et al. 1982), since the time of 
maximum treatment is expected to minimize the propagation effects, and Coulomb 
losses in the solar atmosphere would produce a continuous flattening rather 
than a single break. For small electron events the energy spectra are steeper 
than the spectrum of Figure 3, both below and above the break. Whereas in 
large events, relativistic electrons are occasionally seen above 10 MeV 
(Datlowe 1971), in small events, electrons cannot be seen above a few hundred 
keV. 


Correlations between electrons and protons were studied by Ramaty et al . 
(1980) and more recently by Evenson, Meyer and Yanagita (1981). Figure 4 
(from Ramaty et al. 1980) shows the correlation between 0.5 to 1.1 MeV 
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electrons and 10 MeV protons. As can be seen, for- large events the two 
populations are well correlated, but for smaller events there is an 
overabundance of electrons. This effect is possibly a manifestation of the 
two acceleration phases, with the first-phase producing more electrons than 
protons in small events. The good correlation seen for larger events 
indicates that relativistic electrons just below an MeV and 10 MeV protons 
could be accelerated by the same mechanism. This could be second-phase 
acceleration. The correlation between these two particle population is 
supported by the gamma-ray data. As discussed in I la and seen in Figure 1, 
the >0.35 MeV continuum produced predominantly by <1 MeV electrons, and 
nuclear gamma rays from > 10 MeV protons, have similar time histories. 

Evenson et al . (1981) have recently examined the relationship between 
protons and relativistic electrons at nearly the same energy (' 10 MeV). They 
find that these two particle populations are very poorly correlated. In 
particular, the majority of the proton events have very low (< lO"^) 

electron-to-proton ratios at ' 10 MeV. This is in contrast to the correlation 
seen in Figure 4 where all proton events are accompanied by 0.5 to 1.1 MeV 
electrons. A few of the events, however, show larger electron-to-proton 
ratios, and some of them are as high as 0.2 at 10 MeV. Gamma rays (lines or 
continuum) were seen from all of these events. But it is not clear at the 
present time whether these electron enrichments are caused by the acceleration 
mechanism or whether they reflect different escape condHion for protons and 
electrons. Oppositely directed beams of ions and electrons produced in a 
direct electric field could yield gamma-ray lines at the Sun and relativistic 
electrons in space. Such electric fields are discussed in Section IIIc. 
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Chemical Compositions 

Nuclei heavier than He in solar energetic particles were first detected 
by Fichtel and Guss ( 1961 ) and since then many measurements of such particles 
have been made (see reviews by McDonald, Fichtel and Fisk 1974 , Fan, Gloeckler 
and Hovestadt 1975 , Ramaty et al. 1980 ). While the earlier results indicated 
rough agreement of the energetic particle composition with photospheric 
composition (e.g. Bertsch, Fichtel and Reames 1969 ), with more recent results, 
it became obvious that these two sets of abundances can differ drastically 
from each other. 

The first indications for large abundance anomalies in flare accelerated 
particles came from the observations of Price et al. ( 1971 ) which showed large 
enhancements at low energies of iron-group nuclei over photospheric 
abundances. Further studies (Mogro-Campero and Simpson 1972 , Teegarden, Von 
Rosenvinge and McDonald 1973 ) have revealed such enhancements also for Mg and 
Si. 

The most dramatic departure of a solar energetic particle abundance from 
its photospheric value is that of ^He (Garrard, Stone and Vogt 1973 ). Here 
very large enhancements are occasionally observed in the ^He/^He ratio above 
its likely photospheric value (see Ramaty et al. 1980 for a review of the 
data) . 

Enrichments of ^He in energetic particle populations (for example, the 
galactic cosmic rays) have been generally attributed to nuclear reactions 
between the energetic particles and the ambient medium. But, as first pointed 
out by Garrard, Stone and Vogt ( 1973 ), this interpretation of the solar 3ne 
enrichments, in its simplest form, is inconsistent with much of the ^He 
data. If the 3ne enrichments are due to nuclear reactions of the energetic 
particles, then they should be accompanied by similar enric’iments in ^H and, 
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to a lesser degree, in 3 h. Such enrichments, however, are not observed. 

Several schemes based on nuclear reactions have been proposed to overcome 
this difficulty. These rely on the kinematics and angular distributions of 
the reaction products which favor the preferential escape of ^He (Ramaty and 
Kozlovsky 1974, Rothwell 1976) and the ther-monuclear destruction of and 3 h 
in a model in which the energetic products of the nuclear reactions are 
confined to thin filaments and interact with each other (Colgate, Audouze and 
Fowler 1977). But, as proposed by Fisk (1978) and Kocharov and Kocharov 
(1978), the enhanced ^He abundance in solar energetic particles could te due 
to preferential heating of ambient ^He. Provided that the acceleration 
mechanism has an injection threshold (Section III), such heating would greatly 
enhance the number of accelerated particles. 

Several systematic studies of solar energetic particle composition have 
been carried out recently (McGuire, Von Rosenvinge and McDonald 1979, Cook, 
Stone and Vogt 1980, Mason et al. 1980, Reames and Von Rosenvinge 1981). The 
results of these studies are summarized in Figure 5. Here ratios of energetic 
particle abundances to photospheric abundances are shown for a variety of 
elements as well as for 3ne and ^He. 

For each element or isotope plotted in Figure 5 we consider flares for 
which the energetic particles are rich in 3He and flares from which no 3He is 
observed. The data for these two groups are separated by the dashed vertical 
lines, with 3He-rich flares to the right of the lines and flares with no 3ne 
to their left. 

The closed circles are the data of McGuire, Von Rosenvinge and McDonald 
(1979) who measured the energetic particle composition of eight large non 

■I 

'’He-rich solar flares in the energy range from 6.7 to 15 MeV/nucleon. Such 
compositions were also measured by Cook, Stone and Vogt (1980) from a few to 
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15 MeV/nucleon and by Mason et al. (1980) near 1 MeV/nucleon. The values 
given by crosses and stars are the average particle abundances measured by 
Cook et al. (1980) and Mason et al. (1980), respectively, for their sample of 
large solar flares. As can be seen, these averages are consistent with the 
data of McGuire et al. (1979). 

To the right of the dashed lines are data for 3He-rich flares. The open 
circles represent the recent data of Reames and Von Rosenvinge (1981) for six 
flares, while the values given by the diamonds represent abundances labelled 
by Mason et al. (1980) as anomalous. In particular, diamond-1 is for their 
carbon-poor flares and diamond-2 is for the October 12-13, 1977 flares. Both 
these sets of energetic particles are rich in ^He. 

The distinction between ^He-rich and normal events follows mainly from 
the 3 hc abundance itself. The ^He/O ratio of the former exceeds the upper 
limits on this ratio for the latter by at least an order of magnitude. The 
other abundances, however, are not too different for the two classes of 
events. As can be seen in Figure 5, ^He-rich events can be both rich or poor 
in H (see values of diamonds-1 and 2 for H), although it has been noticed 
(e.g. Ramaty et al. 1980) that ^He-rich flares have on the average low ^H/^He 
ratios. It also follows from Figure 5 that ^He and C tend to be somewhat 
suppressed in 3He-rich events although it is not necessarily true that all ^He 
rich events are C poor. The variability from event to event in ^He, C and 
other elements is larger in 3He-rich events than in normal events. 

Both normal and ^He-rich events show enrichments in heavy elements, 
although these are significantly larger for the ^He-rich events. The 
abundance enhancements in Ca and Fe are in fact almost as large as those of 
^He. The correlation between ^He and Fe enhancements in ^He-rich events was 
first noticed by Anglin, Dietrich and Simpson (1977). 
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Various models have been proposed to account for these abundance 
variations. Meyer (1981) has pointed out that some of the less dramatic 
abundance variation could simply reflect coronal abundances which are 
different from those of the photosphere. Mull an and Levine (1981) suggested 
that differences in Coulomb energy losses (Section Ilia) of partially stripped 
heavy Ions could lower the Injection thresholds of these Ions and hence 
Increase the number of particles that are accelerated. Elchler (1979) has 
pointed out that In acceleration by shocks of finite widths (Section Illb), 
partially stripped Ions with large gyroradll could be preferentially 
accelerated because they are capable of sampling more of the compression at 
the shock. This assumes a rigidity dependent mean free path (Section II Ib) 
and that all Ions have the same Initial velocities. Axford (1982) points out 
that, aside from this effect, the shock may selectively heat Ions of different 
mass-to-charge ratios by plasma processes, so that the Ions will be Injected 
In a manner which depends on this ratio. In Fisk's (1978) model (see also 
Mason et al. 1980), the electrostatic Ion cyclotron waves which heat the ^He 
also heat partially stripped heavy Ions which have A/Z near 3 (e.g. 

56Fe+17) by resonance with the second harmonic of the particle cyclotron 
frequency. Only mechanisms which can selectively preheat certain Ions, such 
as the mechanisms of Fisk (1978) and Kochorov and Kocharov (1978), combined 
with Injection thresholds which are large compared to the thennal speeds of 
the normal Ions In the source, can explain the observed 3He enhancements. 
Isotopic and Ionic Compositions 

In addition to the ^He observations discussed above, the only elements 
whose Isotopic abundances have been observed In solar energetic particles are 
Ne and Mg (Dietrich and Simpson 1979, 1981, Mewaldt et al . 1979, 1981). For 
Ne, It has been found that 22Ne/20j^g _ 3^3 ^ 0.003 and 2lNe/22Ne < 0.1. 
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These abundances are consistent with those of Neon-A In meteorites (Podosek 
1978) which Is believed to represent the primordial Ne Isotopic abundance at 
the time of the formation of the solar system. The Isotopic abundances of Ne 
In the photosphere are not known. In the solar wind, the 22Ne/20ne ratio of 
0.07 ± 0.002 (Geiss 1973) Is significantly lower than In the solar energetic 
particles. The origin of this discrepancy Is not yet understood. But because 
the Mg Isotopic abundances In solar energetic particles, 25Mg/24mg , 
0.14(+0.05, - 0.02), 26Mg/24Mg = 0.15 (+ 0.04, - 0.03) (Mewaldt et al . 1981) 
seem to be consistent with primordial solar system material. It Is possible 
that Ne Isotopic fractionation takes place In the solar wind. Mg Isotopic 
abundances have not yet been measured In the solar wind and are not known In 
the photosphere. 

The first measurements (Gloeckler et al. 1976) of the Ionic states of 
solar energetic particles revealed a charge distribution consistent with that 
of a gas In Ionization equilibrium at about 1 to 2 x 106 k. Subsequent 
measurements (Sclambi et al. 1977) found that the charge states of C and 0 In 
several solar particle events essentially do not vary with energy (from 0.037 
to 1 MeV per charge), In time and from event to event. These results are 
consistent with particles being accelerated In coronal material and traversing 
little material during and after their acceleration. The recent observation 
(Hovestadt et al . 1981) of singly charged He Ions In solar energetic particles 
Indicates that cooler material Is also accelerated. This could be due to 
temperature Inhomogeneities In the corona or the Injection of chromospheric 
material Into the acceleration region. 

Measurements of charge states have Important Implications on the heating 
mechanisms discussed above for producing the ^He and heavy element abundance 
anomalies. Ma Sung et al. (1981) find that 160+5 and 56pe+17 are Indeed 
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present in a ^He rich flare, as would be expected In Fisk's (1970) model. 
However, the existence of appreciable concentrations of such Ions In the 
ambient medium, requires that particles be accelerated from reqlons spanning a 
broad range of temperatures (4xl05 to SxlO^K). 

Recently, Kleckler et al. (1982) have reported that the mean charge of Fe 
during events with large ^He and Fe enrichments Is significantly larger than 
the mean charge of Fe during flares of normal composition. Since the charge 
of Fe during normal events Is not expected to be as high as +17, this result 
Is consistent with the heating model of Fisk (1978) which preferentially 
enhances Fe’’’^^. 


III. MECHANISMS OF SOLAR FLARE PARTICLE ACCELERATION 


Based on the observations described In the previous Section, a solar 
flare acceleration mechanism or combination of mechanisms must fulfill the 
following requirements: the acceleration mechanism must be capable of 
Imparting a large fraction of the available flare energy to energetic 
particles and It must be possible to accelerate a significant fraction of 
ambient particles; the dominant mechanism or the combination of mechanisms 
should be capable of producing the spatial and temporal evolution of the 
different energetic particle populations. In particular the rapid rise of the 
various electromagnetic emissions; the mechanisms must have acceleration 
thresholds, which, with appropriate Injection or preacceleration mechanisms, 
are capable of accounting for the observed compositions and their variations; 
and the mechanisms must be able to produce particles with the observed enprgy 
spectra. 

The observations olso Indicate that the occurence of solar flare Ions In 
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space Is associated with coronal shocks and by Implication with turbulence 
produced by the shock. Therefore, It Is reasonable to Invoke stochastic 
acceleration and direct shock acceleration for solar flare Ion production. On 
an even more basic level, solar flares are associated with the release of 
magnetic energy by the collapse and reconnection of magnetic fields, these 
processes may generate large transient electric fields vhich could accelerate 
Ions and electrons directly. 

In subsection (a) we discuss stochastic acceleration and In subsection 
(b) we discuss shock acceleration with emphasis on recent Ideas of diffusive 
shock acceleration. For both mechanisms we consider energy spectra, the 
problem of Injection, and the acceleration times of the particles. In 
subsection (c) we briefly discuss magnetic reconnection and electric fields 
and their possible role In particle acceleration. 

The results of subsections (a) and (b) are most relevant to Ion 
acceleration In solar flares. The reader is referred to several recent 
reviews of first-phase electron acceleration (e.g. Smith 1979, Ramaty et al. 
1980, Brown and Smith 1980). Our treatment of particle acceleration, as well 
as the treatments that we review, do not discuss the generation and decay of 
the shocks and turbulence, and do not, for the most part, take Into account 
the effect of the accelerated particles on the turbulence or on the shocks 
which accelerate them. This Important aspect of acceleration theory remains a 
subject for future research. 
a. Stochastic Acceleration 

Processes In turbulent plasmas which cause particles to change thelr 
energy In a random way witn many Increases and decreases In energy lead to 
stochastic acceleration. In the original stochastic Fermi mechanism (Fermi 
1949), the process was reflection from randomly moving magnetized clouds. 
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Stochastic acceleration can also result from resonant pitch-angle scattering 
from Alfven waves with wavelength of the order of the particle gyroradius. To 
accelerate paticles these waves must propagate both parallel and anti -parallel 
to the average magnetic field (Skilling 1975). Other modes of stochastic 
acceleration, called magnetic pumping and transit-time damping, occur through 
interaction with magnetosonic waves whose wavelengths are much longer than the 
particle gyroradius (Kulsrud and Ferrari 1971, Melrose 1980, Achterberg 
1981). These modes require additional pitch-angle scattering to keep the 
particles isotropic. Langmuir (plasma) waves or other electrostatic waves 
with phase velocities of the order of the particle speed will also accelerate 
particles stochastically (Melrose 1980). 

When the random energy increments are small compared to the particle 
energy, stochastic acceleration results in a diffusive current in momentum 
space, Sp = - Dpp3f/3p, where y- is the magnitude of the momentum, f(p) is the 
number of particles per unit volume in phase space and Sp is measured in cm”'^ 
momentum “2sec-l. Particles injected at some momentum Pp win diffuse in 
momentum to larger and smaller p. In terms of f, the differential particle 
intensity per unit energy per nucleon is given by dJ/dE = A p2f, where A is 
the nuclear mass number. Additional non-diffusive energy changes can be added 
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(3) 


where dp/dt represents convection in momentum space due to processes which 
change the energy of all particles (e.g. ionization or Coulomb losses). 

The momentum diffusion coefficient Dpp depends on the nature of the 
stochastic process. If the process is hard sphere scattering with mean free 
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path X, then Dpp can be derived (Parker and Tidman 1958, M. A. Lee, private 
communication 1978) from the Boltzmann equation. This yields 


Dpp = p^6V)^/3vX 


(4) 


where (6V)2 is the mean square velocity of the scatterers and v the particle 
speed. 

If the stochastic acceleration is due to resonant pitch-angle scattering 
from Alfven waves, the momentum diffusion coefficient obtained from quasi- 
linear theory is (Skilling 1975) 


2 2 
2P Va 
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(5) 


where the Alfven speed, D+(D_) is the pitch-angle scattering coefficient 

due to forward (backward) propagating Alfven waves, and p is the cosine of the 
particle pitch angle in the mean field. It is clear from this equation that 
there is no stochastic acceleration due to Alfven waves unless the waves 
propagate in both directions. This requirement occurs because the electric 
fields of Alfven waves in one direction can be Lorentz- transformed away and so 
cannot accelerate particles. For example, Alfven waves generated by the 
streaming of energetic particles (Wentzel 1969) propagate only in the 
direction of the streaming and hence do not accelerate the particles. 

By comparing equations (4) and (5) we can define an effective mean free 

A ? ? A 

path for stochastic acceleration by Alfven waves, x = p Vw3vD . This 
xA is, in general, a function of particle rigidity which is determined by the 
power spectrum of the Alfven waves. Let W^(k) be the energy density per wave 
number k in waves propagating in the + or - directions. Then, from quasi- 
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linear theory (Hasselmann and Wibberenz 1968, Jokipii 1971, Luhtnann 1976), 


D^(ii) = v(2itZe/pc)^ W^(l/ur^){l-u^)/lul, 


( 6 ) 


where Ze is the particle charge and its gyroradius in the average magnetic 
field, B. For example, if 0+ = 0_ and W{|() “ kn-2^ 


X^(R) = ^ 


(2-n) (4-n) « R^ 


96Tt"W(l/r^) 


(7) 


For stochastic acceleration due to long-wavelength magnetosonic waves 

(assuming adequate particle pitch-angle scattering) the expression For Dpp 
from quasi -linear theory is 


D 



= i 




<k> 


( 8 ) 


p 

(adapted from Achterberg 1981), where <6B > is the mean square of the 
fluctuations in the field magnitude and <k> is the mean wavenumber of the 
magnetosonic waves. The number ^ depends on the angular distribution of the 
waves and it is usually assumed that ^ ~ 1. Note that the corresponding 
to equation (7) is B <k>” /<6B > and is independent of particle momentum or 
charge. It is not necessary for the magnetosonic waves to propagate in both 
directions to accelerate particles, but as already mentioned a certain level 
of pitch-angle scattering is required to isotropize the particles. This 
condition is D^+0 (Achterberg 1981). Unlike for acceleration by 

Alfven waves, the waves which do this scattering are not required to propagate 
in both the + and - directions. 
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The relative importance of acceleration by magnetosonic and by Alfven 

waves is given by D In order of magnitude, this ratio is 

PP PP 



where k = r is the resonant wave number and is the total energy 
r c 

density in the long wavelength magnetosonic waves. Since by assumption 

<k>r <<1, acceleration by Alfven waves dominates when their energy density is 
c 

comparable to that in magnetosonic waves, but only when there is appreciable 

power in Alfven waves propagating in both directions. 

The momentum diffusion coefficient for isotropic Langmuir 

turbulence, D is given by Melrose (1980; e.q. 8.13). This results in 
op 

L 2 

a X « (A/Z) times a function of particle velocity. 

The acceleration models which we now discuss in detail consider the 
physical mode of stochastic acceleration only through the momentum diffusion 
coefficient Dpp. Neglecting spatial convection, averaging over some volume of 
space and introducting an escape time T from this volume, particle 
conservation results in the transport equation 

P 

where Q(p,t) is the particle source in momentum space. Equation (10) has been 
applied to the acceleration of solar energetic particles by Barbosa (1979), 
Ramaty and Lee (Ramaty 1979) and Mull an (1980). 

In Barbosa's (1979) model (hereafter referred to as model B) the 
acceleration is by Alfven waves with a power law spectrum for Vl(k), T is 
assumed to be a power law in velocity, T«v''’ , and dp/dt is neglected. In 
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Mullan's model (1980, referred to as model M), the scattering Is from neutral 
sheets in turbulent motion leading to a constant X; numerical integration is 
used in order to incorporate realistic losses due to Coulomb scattering and 
adiabatic deceleration in an expanding turbulent region, followed by an 
energy-independent escape. The Ramaty-Lee model (Ramaty 1979, model RL) is a 
special case of the other two, since it assumes that both the mean free path 
and escape time are energy-independent and that dp/dt is negligible when 
acceleration is effective. These are rather special assumptions which need to 
be physically justified. Nevertheless, the predicted energy spectrum of this 
model fits the interplanetary observations of proton and a-particle spectra a1 
lAU very well (see Section Ilb). Because of this and because of its inherent 
simplicity, we proceed to first discuss the RL model. 

With a steady source of q particles/cm^sec at momentum Pq, equation (10) 
becomes 




q6(p-p^) 



> 


( 11 ) 


2 

where the acceleration parameter a = (5V) /xc. The steady-state solutions 
of this equation are characterized by the dimensionless, energy-independent 
constant aT. For ul trarelativistic particles 


P -3/2 ± (9/4 + 3/aT) ^2 

f = ^ , 

47rp^\ (1 + 4/3aT) ^2 




( 12 ) 
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where the plus sign applies to p<po and the minus sign to p>Po* 
non-relativistic particles and p>p^j 


f = ^3—2 ^ Ip (2(3p /mcaT) )K (2(3p/mc»T) ) , (13) 

4irp^ mco ^ 

where m is the ion mass, and I 2 and K 2 are the modified Bessel functions of 
order 2 (Abramowitz and Stegun 1966). The arguments of I 2 and K 2 ®re 
interchanged for p<Pq. The corresponding intensity per energy per nucleon 
is dJ/dE = Ap2f. 

The spectrum of dJ/dE is not a power-law. It has an energy-dependent 
slope y(E) = - dln(dJ/dE)/dlnE which approaches zero at low energies. From 
the asymptotic expression for K 2 (x) « x "^'^exp(- x), we find that at energies 
E>> 3.26 a2f2 MeV/nucleon 


dJ 

W 


oc 


E^/®exp {- (E/(3.26(aT)^))^/M . 


(14) 


This spectrum steepens with increasing energy up to the fully relativistic 
domain where equation (12) applies and y(E) = (l/2)(9 + 12/aT)^'^- 1/2. 

As discussed in Section lib, McGuire, et al. (1981) have found that 
equations (13) or (14) fit the spectra of protons and a-particl<s in 
interplanetary space quite well with aT around 0.025 for protons. Different 
values of aT, however, are required for these two particle species; in many 
events the values of aT which fit the helium spectra are less than those for 
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protons. The solar energetic particle spectra, therefore are not strictly 
velocity -depen dent as Implied by a constant x. The possibility of a rigidity 
dependent diffusion mean free path Is considered In the model of Barbosa 
(1979) which we discuss next. 

As already mentioned, the B model uses the quasi -linear theory of 

particle resonant scattering to evaluate Dpp ftom a given spectrum of Alfven 
waves. The steady state spectrum of non-relativistic particles In this model, 
with x«*R'^(cf. equation 7) and T«v? Is 

^ ( 3 ^, ^ H5) 

^A ^ 


where s = (n+l-m)/2>0 and v = (2-n)/2s. By using the asymptotic expression 
for K^{x), equation (15) can be written as 


dJ _(3+n+m)/8 


( /A.n/2 
exp (-( 7 .) 


(- 


E >(l+n-m)/4. 

2J i 


3.26K 


( 16 ) 


which 1 s valid if the argument of the exponential Is much larger than unity. 
The constant K In equation (16) for the B model Is equivalent to aT In 
equation (14) for the RL model. For n = m = 0 equations (15) and (16) reduce 
to equations (13) and (14), respectively, I.e. the B model reduces to the RL 
model . 

For n>0 equation (16) predicts a-particle spectra, which, when expressed 
as functions of energy per nucleon, are steeper than the proton spectra for 
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the same K. This feature of the B model, not present In the RL model, 1s 
qualitatively consistent with the observations, as discussed above. But the 
comparison of equation (16) with the solar flare proton spectra for a variety 
of m's and n's shows that the best fit Is obtained for m = n = 0. Other 
values make the slope of the spectrum vary more rapidly with energy than 
observed. For example, If m = 0, we find that equation (16) fits the observed 
proton spectrum shown In Figure 2 only if n ^ 0.1. For such small values of 
n, the a-particle spectrum is only slightly steeper than the proton spectrum 
with the same K, while the observed proton and a-part1cle spectra shown In 
Figure 2 are considerably different. Thus, while the B model allows rigidity 
dependent spectra, the fit of the calculations to the data Implies that this 
rigidity dependence is not very pronounced. 

In the M model, the scattering elements are assumed to move with 
velocities 5V of the order of shock velocities In the corona. Since x Is 
assumed to be energy independent Dpp, 1$ given by the same expression as In 
the RL model, equation (4). The dp/dt term in equation (3), however. Is taken 
Into account. In particular a time dependent combination of Coulomb and 
adiabatic deceleration losses followed by energy Independent escape Is used. 
Because of this complexity, the resultant energy spectra can only be given 
graphically. They are approximately constant below 1 MeV /nucleon, very steep 
above 100 MeV/nucleon, and have y ^ 3 at ^ 50 MeV/nucleon. This Is In general 
agreement with observed spectra, although detailed comparisons over an 
extended energy range have not yet been made. 

A very Important question In all particle acceleration theories. 

Including stochastic acceleration. Is that of Injection. We first note that 
the basic concept of stochastic acceleration assumes that the energy changes 
are small compared with the particle energy and therefore the particle 
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velocity must be much greater than 6V. Furthermore, for resonant scattering, 

ions must have v>V^ to scatter from Alfven waves and electrons must have 
v>43V;^ to scatter from whistlers (Melrose 1974). 

An additional injection condition is set by the requirement that the 
systematic acceleration rate due to diffusion in momentum space be larger than 
the ionization and Coulomb energy loss rates of the particles. The systematic 
acceleration rate is (e.g. Ramaty 1979) 




(17) 


and for Dpp from equation (4) 



+ 2Mc2 ) ) 


(18) 


Here M and E are the proton mass and kinetic energy per nucleon for nuclei, 
and the electron mass and kinetic energy for electrons. 

Energy loss rates due to ionization in a neutral medium and Coulomb 
losses in a fully ionized medium were summarized by Ramaty (1979). These loss 
rates together with the systematic acceleration rates for nuclei and electrons 
(equation 18) are plotted in Figures 6 and 7 for neutral and ionized media, 
respectively. Particles can be accelerated only if the rate of systematic 
energy gain exceeds the rate of energy loss. Depending on the ratio of the 
acceleration efficiency, a, to the ambient density, n, the energy gain curve 
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maiy or may not Intersect the energy loss curve. In the former case an 
Injection mechanism is required which preaccelerates the particles to E - E^, 
where Eq energy where the two curves Intersect. 

The values of a/n Indicated In Figures 6 and 7 were chosen such that E^ 
for electrons Is around 0.1 MeV, an energy at which there seems to be a 
transition from first phase to second phase acceleration (Section II). For 
protons, however, these a/n's are such that the systematic gain Is larger that 
the loss at all energies, and in this case It Is possible to accelerate 
ambient particles directly. However, stochastic acceleration still requires 
that the particle velocities be larger than 6V. Assuming that 6V Is of the 
order of the Alfven velocity, the particles must be accelerated either from a 
high B plasma (B = particle pressure/magnetic pressure) or possibly from the 
non-thermal tall of a low-B plasma. Such non-themal tails are observed In 
the solar wind (Ogllvie, Scudder and Olbert 1978, Scudder and Ogllvie 1979). 

If the Alfven speed Is the threshold, then 3ne-r1ch events cannot originate In 
h1gh-B plasmas. 

The solar energetic particle abundances and their variability must be 
strongly dependent on the Injection process as we have discussed In Section 
Ilb. Injection questions have not yet been fully Investigated and they remain 
outstanding problems for future research. 

The acceleration time of stochastic acceleration can be studied from the 
time dependent solutions of equation (10). Such solutions were obtained In 
the RL model for impulsive Injection of particles at t » 0 and p = p^^ 
with T » and \ = const. In Figure 8 (from Ramaty 1979) we show the 
differential proton and electron number densities (N(E) = 4itp^f/v) at various 
times t after Injection. As can be seen, In ? given time, protons are 
accelerated to much higher energies than electrons. This result Is the direct 
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consequence of the acceleration rates shown In Figures 6 and 7 which assume 
the same a for protons and electrons. 

From the gamma- ray observations (Section Ila) It follows that In at least 
some flares protons and nuclei are accelerated to energies greater than 10 
MeV/nucleon In about 1 sec. As can be seen from Figure 8, acceleration of a 
significant fraction of the protons to suqh energies requires 
that <xx - 0.1. Then If t » 1 sec, a * 0.1 sec"^. Because 6V can probably 
not be greater than 3x1 08 cm/sec, X must be less than 3xl07cm. Since this 
value 1s 50 times the gyroradlus of a 20 MeV proton In a 1 Gauss field, we 
could expect such short mean free paths In the strong turbulence of solar 
flares. But there is no direct observational evidence for such turbulence. 

The characteristic acceleration time in the B model Is of the same form 
as In the RL model. Barbosa (1979) estimates that the acceleration time to 
10 MeV can be as short as 10 seconds. As above, this requires strong 
turbulence on scales lO^cm. Mull an (1980), on the other hand, estimates 
from numerical solutions that the time to accelerate keV protons to '• 5 MeV Is 
about an hour. This time is too long for the acceleration of the gamma-ray 
producing protons and nuclei, but It could be adequate for the acceleration of 
the solar flare particles observed in interplanetary space. 
b) Shock Acceleration 

Solar flare shocks propagate upward through the solar corona at speeds of 
about 500 to 2000 km/sec, as indicated by Type II radio bursts (Smith and 
Goldman 1982 and Section Ila), and laterally through the chromosphere where 
they are seen as Moreton waves (Uchida 1974). The occurence of solar- 
energetic particles in space is strongly correlated with flares having Type II 
bursts (Svestka and Fritzova 1974). Flare shocks are observed to accelerate 
particles in interplanetary space (Richter and Keppler 1977), as are 
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co-rotat1ng shocks (Barn?s and Simpson 1976, McDonald et a1. 1976) and 
planetary bow shocks (Asbrldge, Bame and Strong 1968, Lin, Meng and Anderson 
1974, Scholer et a1. 1980, Zwicki et a1. 1981). A flare shock can transport 
particles In an energy-independent manner through the corona until they escape 
onto open field lines. Shock acceleration has been recently reviewed by 
Toptygin (1980) and Axford (1982) and applied to solar flares by Achterberg 
and Norman (1980), Decker, Pesses and Armstrong (1981) and Lee and Fisk 
(1982). See also McDonald (1981) for review. 

There are basically two types of shock acceleration: scatter- free. In 
which particles gain energy by reflection 1n a single shock encounter 
(Sonnerup 1973, Chen and Armstrong 1975, Pesses, Decker and Armstrong 1982) 
and diffusive, 1n which particles gain energy by repeated scattering between 
the converging plasmas on either side of the shock (e.g. Axford, Leer and 
Skadron 1977, Axford 1982). The scatter- free mechanism can enhance the 
particle energy by about an order of magnitude If the shock Is nearly 
perpendicular (I.e. the magnetic field Is nearly perpendicular to the shock 
normal), but In that case, only particles with speeds which are already much 
greater than the shock speed can be reflected. Further acceleration, however, 
requires multiple reflections. These are possible If there Is particle 
scattering In the fluid flow or If the particles are trapped between 
converging shocks In a flare loop (Wentzel 1965). 

Acceleration by diffusive scattering across the shock Is a first order 
Fermi process. In the sense that every shock crossing results In an energy 
gain. It Is In principle more efficient than stochastic acceleration because 
It derives energy directly from the compression of the flow at the shock. For 
this mechanism to be effective, however, there must be adequate particle 
scattering both upstream and downstream of the shock. The passage of the 
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shock Is expected to generate turbulence downstream which will scatter the 
particles. Scattering upstream* however*1s more problematic (Holman, lonson 
and Scott 1979). Observations (Tsurutani and Rodriquez 1981) of 
Interplantetary shocks and planetary bow shocks show that when they are nearly 
parallel there Is a very turbulent foreshrick region capable of scattering 
particles. Such a region could be produced by the accelerated particles 
themselves (Achterberg and Norman 1980, Lee 1982). 

In the simplest example of a plane shock where the only losses are due to 
convection of the particles away from the shock downstream, the energetic 
particle density In phase space Is given by a power law, f « p-3VMV^ where V 
Is the shock speed and aV the discontinuity In the plasma speed at the 
shock. In terms of the compression ratio, r, at the shock, V/aV * r/(r-l). 

For a strong shock In a non-relativistic fluid, r » 4 and hence f(p) « p"^. 

For weaker shocks, 4>r>l and tk -efore the power-law Is steeper. Blandford 
and Ostrlker (1980) have modeled cosmic-ray acceleration and propagation In 
the Galax/ on the basis of this result. Including acceleration by shocks with 
a spectrum of compression ratios. They have been able to reproduce the 
observed cosmic- ray spectrum which approximates a power law over a wide range 
of energies. 

Unlike galactic cosmic rays, solar flare particles do not have power-law 
spectra (e.g. Figure 2). In fact, none of the energetic particle populations 
which are observed to be accelerated by shocks In Interplanetary space have 
power-law spectra. Diffusive shock acceleration produces such curved spectra 
rather than power laws when particle energy losses or escape losses are 
significant, or when the shock has a finite size or lifetime, compared to the 
natural scales of the shock acceleration. For example, diffusive escape 
losses from the finite bow shock of the Earth (Elchler 1981, Forman 1981b, 



40 


ElUson 1981a, Lee 1982) and adiabatic losses near interplanetary shocks (Fisk 
and Lee 1980, Forman 1981a) have been shown to reproduce the observed spectra 
quite well. 

In the present subsection we describe the general methods used to obtain 
energetic particle spectra from shock acceleration, and provide a simple 
solution of the appropriate transport equation which provides a spectrum that 
fits the solar flare data quite well. 

The transport equation which describes diffusive shock acceleration is 
similar to equation (10) without the stochastic acceleration term proportional 

to Dpp^ but with additional terms due to convective transport, spatial 
diffusion and adiabatic compression of particles in the plasma flow. This 
equation is given by (Axford 1982) 
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Q(p,r,t), (19) 


where f, T and dp/dt have been defined in connection with equation (10), ¥ is 
the plasma velocity and k the spatial diffusion tensor which couples the 
energetic particles to the plasma converging at the shock. The terms 
containing V and k are essential for the description of acceleration across 
the shock front. The Injected particles are explicitly introduced in equation 
(19) by the source term Q; the injection may also be treated as a boundary 
condition such that f approaches a given value fg far upstream. The losses 
due to particle escape can be treated via the escape time T, or as diffusive 
escape. In the latter case the scattering becomes negligible at a finite 
distance from the shock. 

The usual method for deriving a steady-state » 0) particle spectrum 
is to first solve equation (19) separately on each side of the shock and then 
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to match the two solutions at the shock by imposing boundary conditions. These 

conditions are that both the energetic particle density and the normal 

2 af 

component of the spatial streaming ( S = - 4np (Vp -^ + k»vf)) of these 
particles be continuous at the shock. Toptygin (1980) has shown that this is 
an appropriate approach even though equation (19) is not valid very close to 
the shock. 

We present here a steady-state solution of equation (19) by assuming that 
dp/dt is negligible and T is constant (as in the RL model for stochastic 
acceleration), that only diffusion along the shock normal is impor*tant (hence 
k can be treated as a scalar), that the shock is an infinite plane at x = 0 
and infinitely thin, and that the particles are steadily injected at the 

7 7 

shock, i.e. Q = q 6(x)5(p-p^)/(4TTp^ ) where q is measured in particles cm 

sec For these conditions, the solution on either side of the shock is 

given by 


f(x,p) = f(0,p) exp (-B^lxl) , 


or POOR 


( 20 ) 


where upstream 

= (V + (V^ + 4k^/T)^^)/2k^ (21) 


and downstream 


= B 2 = (-(V - AV) + ((V - AV? + 4k2/T)^'^)/2k 


2 ‘ 


( 22 ) 


The streaming continuity at the shock gives the equation 
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(23) 


The exact non-relativistic solution of this equation Is 


f(0,p) (-3(VI, + (V - AV)Ip)/AV] . (24) 

4irp ^0 
0 

where 

I. = n+n.)^^2_i _ ctnh'^ (1+n^)^^- In (n^/4)^'^ , (25) 

with n^- = 4x^. p/(3mV?T) . Since for n^-+0, the distribution function Is a 
power law for v << 3V^T/4x. with the same spectral index (3V/AV) as for the 
case of a plane shock with no losses. However, when v becomes comparable to 
the lower value of 3V^T/4x^. , the shape of the spectrum Is determined by the 
exponential term In equation (24). In other words, the spectrum Is determined 
by the losses on the side of the shock where V^/x^. Is lower. This Is probably 
the upstream side. 

Assuming for simplicity that V^/x is the same on both sides of the shock, 

and that 3 V/aV = 4, as for a strong shock, the differential Intensity per 

2 

energy per nucleon, for v>> 3V T/4x, Is 

exp (-(E/(0.422(aT)^))^/4), (26) 

where aT = V^T/(xc). Equation (26) for shock acceleration Is quite similar 
to equation (10) for stochastic acceleration In the RL model. These two 
spectra are compared In Figure 9 using values of aT which fit typical data In 
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both cases. As can be seen, the difference between then is most pronounced at 
the lowest energies where the solar flare spectrum is unknown because of 
adiabatic energy losses in the interplanetary medium (Section IV). 

As for stochastic acceleration, for shock acceleration as well, the 
question of injection is very important. Ionization and Coulomb energy losses 
to the ambient medium have the same role in determining injection conditions 
in shock acceleration as they do in stochastic acceleration. In addition, for 
diffusive shock acceleration, particles downstream must have sufficient 
velocity to overtake the shock. This is at least (V - aV) directed towards 
the shock, and increases as cos-J^ij;, where ii; is the angle between the 
downstream field and the shock normal. The velocity V-aV is at least as great 
as V/y 0^' 3nd with the additional cos"2\|) factor, the threshold for shock 
acceleration is expected to he higher than for stochastic acceleration. 

Another injection condition is set by the finite width of the shock which 
could depend on many parameters including the pressure of the accelerated 
particles. When this pressure is taken into account (Axford, Leer, and 
Skadron 1977, Eichler 1979, Drury and Volk 1981, Drury, Axford and Summers 
1932) all or part of the velocity change AV is smoothed out over a length 
scale '• R/v, where R is the diffusion coefficient of the energetic particles 
averaged over their energy spectrum. As pointed out in Section lib, this 
smoothing is expected to affect the composition of the accelerated particles 
(Eichler 1979). Drury, Axford and Summers (1982) show analytically that when 
k is independent of energy, the smoothing causes the dominant accelerated 
species (i.e. protons) to have a steeper spectrum than in the case of an 
infinitely thin shock. Minor species which are only partially stripped could 
have larger diffusion coefficients than the protons if the diffusion mean free 
path is rigidity dependent, and therefore for them k>R. Drury, Axford and 
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Summers (1982) show that the spectrum o', such minor species is flatter than 
for protons and approaches that of an ■‘nfimtely thin shock. 

Ellison (1981b) has studied these effects with a Monte Carlo calculation 
appropriate for a parallel shock. In this treatment, the plasma flow and the 
energetic protons are required to conserve mass, momentum and energy flux, the 
mean free path is assumed to be proportional to gyroradius, and particles are 
injected from the shocked plasma downstream in which the temperatures of 
different ions are proportional to t'.eir mass. Ellison (1981b) finds that the 
shock is indeed broadened and t.iiat inere are modest enhancements of energetic 
particles which increase with the mass-to-charge ratio. 

The acceleration time shock acceleration can be obtained from the 
time-dependent solution of equation (19). Such solutions have been obtained 
for various initial vonditions and geometries by Fisk (1971), For-man and 
Morfill (1979) and Toptygin (1980) and have been reviewed by Axford (1982). 

The general result is that 


dE _ aV r\ . '‘2> 



(27) 


This expression shows explicitly that if the mean free path on either side of 
the shock is very large, the rate of energy gain is very small. Since, as 
already mentioned, shocked gas downstream is expected to be turbulent, the 
efficiency of of the acceleration depends critically on the presence of 
turbulence upstream. 

In the absence of upstream turbulence, particles can still be accelerated 
stochastically by the shock-generated turbulence downstream. By comparing 
equations (18) and (27) we find that the ratio of the shock acceleration rates 
and the stochastic acceleration rate downstream is of the order 
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VaV/(6V)2 + x^)), where 6V is the velocity of the turbulent elements 

and X2 and x^ are the downstream and upstream mean free paths, respectively. 
Thus, if the mean free paths are comparable, shock acceleration will be more 
rapid since V and aV are expected to be larger than 6V. 
c Acceleration in Direct Electric Fields 

The previous two sub-sections have dealt with mechanisms that can be best 
applied to the acceleration of ions in solar flares. In addition to 
stochastic and shock acceleration, there is also the possibility of 
accelerating particles in direct electric fields. Such fields are associated 
with magnetic reconnection in the vicinity of magnetic neutral points and 
current sheets (Syrovatskii 1981), where they appear perpendicular to the 
magnetic field. Particle acceleration is also possible in electric fields 
parallel to the magnetic field (Colgate 1978). Parallel electric fields arise 
from the interruption (due to plasma instabilities) of the parallel currents 
associated with twisted magnetic flux tubes and from the formation of double 
layers of electric charge (e.g. Spicer 1982). 

The application of direct electric fields to the acceleration of 
nonrelati vistic electrons in solar flares has been discussed most recently by 
Spicer (1982). (See also Smith 1979, Heyvaerts 1981). Here we wish to 
emphasize the role that direct electric fields could play in the very rapid 
production of ions and the acceleration of relativistic electrons. 

As discussed in subsections (a) and (b), both stochastic and shock 
acceleration could accelerate ions quickly enough to account for the gamma-ray 
observations, but it is not clear that the mean free paths are, in fact, 
sufficiently short and the turbulent and shock velocities are large enough to 
account for the rapid acceleration. On the other hand, in an electric field 
model (such as that of Colgate 1978), the acceleration time could be as short 



46 


as 0.1 sec for a loop of length *>• 109cm. 

Relativistic electrons around 10 MeV, however, cannot be produced fast 
enough and in sufficient quantities in stochastic and shock acceleration 
because the particle energy gain per collision is proportional to particle 
mass. But there seems to be a good correlation between interplanetary 
electrons of such energies and ganwna-ray flares (Evenson, Meyer and Yanagita 
1981; Section Ilb). Since the gamma rays are produced by ions of 19 
MeV/nucleon, this correlation could indicate electron and ion acceleration to 
the same energy by an electric field. 

In this Section, we review particle acceleration processes in current 
sheets and in the vicinity of magnetic neutral points. For discussions of 
acceleration by parallel 'lectric fields, the reader should consult Colgate 
(1978), Smith (1979), Heyvaerts (1981), and the review of Spicer (1982). 

Current sheets and magnetic neutral points are involved in a wide variety 
of contexts in solar physics. The primary energy release in solar flares may 
occur by magnetic reconnection (e.g. Kahler et. al . 1980, Sturrock 1980). 
Parker (1979 and references therein) has emphasized the near-singular 
conditions required for magnetostatic equilibrium, and has proposed that the 
reconnection of tangled magnetic fields occurs routinely and inevitably in 
nature. Syrovatskii (1981) has recently reviewed the formation of current 
sheets. In view of the important role suggested for current sheets and 
reconnection in the solar atmosphere it is natural to consider particle 
acceleration mechanisms associated with these phenomena. Reconnection and 
associated particle acceleration have been extensively studied in the Earth's 
magnetosphere (Stern and Ness 1982). 

In models of steady state reconnection, magnetic fields of opposite 
polarity are brought together by a flow of fluid. A detailed discussion of 
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the fluid theory of steady state reconnection has been given by Vasyllunas 
(1975). Early descriptions of the reconnection process and Its application to 
solar flares were given by Parker (1957, 1963) and Sweet (1958). Two of the 
best known examples of reconnection are the models of Petschek (1964) and 
Sonnerup (1970) which Involve x-type neutral points. The magnetic field lines 
are sketched In Figures 10a and 10b for these two models. These field lines 
lie in a plane (say, the x-y plane) and there Is a constant electric field, E, 
perpendicular to that plane. The electric field is related to the magnetic 
field B, flow velocity V, and current J In the plasma by Ohm's law 

E + V X B/c *- nJ. (28) 

where n is the resistivity. The Inclusion of other terms In equation (28) Is 
discussed by Vasyllunas (1975). The resistive term 1s unimportant everywhere 
except in a small region near the x-type neutral point where B becomes small 
and J Is large (because the gradients of B are large). Here the frozen In 
condition of the magnetic field Is violated. Rather than being transported by 
the fluid, the field diffuses through the semi-stagnant plasma. The region, 
therefore. Is referred to as the diffusion region. 

As shown In Figure 10, cold fluid flows Inward along the x axis. Most of 
the fluid never enters the diffusion region, but passes through a shock front 
(or pair of fronts. In Sonnerup 's model) and flows out parallel to the y 
axis. In Petschek 's model, the Inflow occurs at a small fraction (~«1) of the 
Alfven speed, while the outflow Is exactly Alfvenic. This acceleration takes 
place at the expense of magnetic energy. In Sonnerup's model, the inflow can 
be faster than In Petschek 's model. In which case, less acceleration of the 
bulk flow occurs. In both models, however, the thermal speed of the Ions 
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increases across the shock front to about the Alfven speed. The increase of 
electron thermal energy depends on the structure of the shock and the nature 
of the coupling between electrons and ions. 

Similar results have been obtained by Hayashi and Sato (1478) who have 
simulated time dependent, compressible reconnection by solving the fluid 
equations numerically. In their model, an initially planar neutral sheet is 
compressed into an x-type neutral line by imposed fluid flow directed inward, 
perpendicular to the initial field lines. They assume that when the current 
density J exceeds some critical value Jc» anomalous resistivity develops and 
increases as a power of J - Thus, reconnection of the magnetic field- 
lines and Joule beating of the plasma near the neutral line begin when the 
neutral sheet has been sufficiently compressed by the incoming flow. They 
find that, as in the analytic models of reconnection, the regions of inflow 
and outflow are separated by shocks and that the outflow speed is on the order 
of the Alfven speed. 

If ions gain flow velocities on the order of the Alfven velocity, their 
energy per nucleon is 2.5xl04 B2/n, where B is in gauss and n in cm-3. While 
for most combinations of B and n in the solar atmosphere this energy is not as 
high as the observed particle energies, this bulk flow acceleration could 
nonetheless be important as an injection mechanism. In certain exceptional 
cases, however, B could be high enough and n low enough for the particles to 
achieve energies as high as 10 MeV/nucleon (e.g. Sonnerup 1973) which would 
suffice, for example, for the production of nuclear gamma rays. 

In addition to the acceleration and heating of the plasma as it passes 
through the reconnection region, the possibility exists of direct particle 
acceleration in the electric field at the magnetic null line. 

An example of such acceleration was discussed by Speiser (1965). He 
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considered a magnetic field with field lines parallel to the y axis with field 
rever sal at x = 0, 
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where ty is the unit vector in the y direction and Bj is a constant. Here the 
field equals Bj for |x|>d and goes linearly to zero for |x|<d. A constant 
electric field E = VB.e^/c is also present and assumed to be continuous at 
x=0. Everywhere in space except near the field reversal region at x=0, the 
particle orbits consist of gyration in the x-7 plane (the usual Lar-mor motion) 
combined with the ExB drift in the x direction, toward the plane x=0. Since 
for lxl>d the motion in z is oscillatory, the average energy gain is zero in 
this region. However, once a particle has drifted to near the x=0 plane where 
the field is very weak, it does not gyrate - it is simply accelerated in the z 
direction. The final energy, E, of the particle depends on the extent, L, of 
the current sheet in the z direction. Numerically, E in MeV equals 10"^ BLV, 
where B is in gauss, L in km and V km/sec. Thus particles can be accelerated 
to tens of MeV if B,V and L are large enough. 

Whether neutral sheets have sufficient lengths, and, whether or not 
particles can remain in the sheet for such lengths, are problems which have 
not yet received clear-cut answers for solar flares. 

In order to estimate the distance the particles stay in the neutral 
sheet, particle orbits in various magnetic geometries have been calculated by 
Speiser (1965), Friedman (1969), Bulanov and Sasorov (1976), Bulanov (1980), 
and Syrovatskii (1981). The results depend strongly on assumed geometries, 
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and analytic treatjnents can only be made for the simplest field 
configurations. Time dependent electric and magnetic fields have also been 
considered (Burke and Layzer 1969, Levine 1974). Levine extended Burke and 
Layzer's work by trying to estimate the effect ov Coulomb collisions between 
the test particles and field particles In the ambient medium. The existence 
of an energy loss mechanism which competes with acceleration sets an upper- 
limit on the time constant for collapse of the neutral sheet such that 
particle acceleration can still occur. Recently, Mullan and Levine (1981) 
have studied the Implications of the collapsing magnetic neutral sheet model 
for the composition of solar flare accelerated Ions (see Section Ilb). 

Attempts have also been made to calculate the energy spectrum of 
particles accelerated In neutral sheets. In the absence of stochastic 
processes such as Coulomb collisions or wave-particle Interactions, the phase 
space distribution function f of the test particles satisfies Llouville's 
theorem; the density of particles In phase space Is constant on phase space 
trajectories. Therefore, given an Initial spatial distribution of Injected 
particles and solutions for the particle orbits, the momentum distribution 
function of particles leaving the acceleration region can be calculated. 
(Bulanov and Sasorov 1976, Bulanov 1980). They find that when the Initial 
distribution of particles Is uniform In space, the spectrum Is an exponential 
In energy or In a fractional power of the energy. Energy spectra have also 
been calculated numerically by Bulanov (1980) who finds good agreement with 
the analytic results. Friedman (1969) has used numerical techniques to 
calculate final energies for test particles In Petschek's (1964) reconnection 
model. A generally recognized problem with particle acceleration In neutral 
sheets Is that only very few particles are accelerated. 
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IV. SOLAR FLARE ENERGY SPECTRA IN INTERPLANETARY SPACE 


The spectra of particles observed at any one time In Interplanetary space 
are not exactly the same as the spectra released at the Sun because of 
velocity dispersion caused by energy-dependent diffusion, and because the 
solar wind convects and decelerates energetic particles as it expands. If 
these solar wind effects are small, the effect of velocity dispersion can be 
eliminated in prompt, magnetically well-connected events by using the time-of- 
maximum (TOM) energy spectrum forned by the maximum intensity at each eneroy 
(Lin 1970, Van Hollebeke et al . 1975). We justify this method below. 

At energies less than about 1 MeV/nucleon for ions and 1 keV for 
electrons, convection and adiabatic deceleration will make the TOM spectrum 
flatter than the spectrum at the Sun, but the precise effect on the spectrum 
is uncertain. In addition, the TOM spectrum becomes less reliable for such 
low-energy particles because their propagation becomes more sensitive to 
inhomogeneties in the solar wind and magnetic fields (e.g. shocks). 

The important questions which need to be resolved before interplanetary 
spectra can be used to test models of the acceleration and release process at 
the Sun are: What is the most convenient way to estimate the spectrum at the 

Sun from spacecraft observations of energetic flare particles; above what 
energy is that spectrum reliable; and what is the approximate nature of the 
change in the spectrum below that energy. First, we show why the TOM spectrum 
in simple well-connected events is a good approximation to the spectrum at the 
Sun. 

The observed f(p,r,t) is related to the spectrum injected into 
interplanetary space by the transport equation. 
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(30) 


which Is analogous to equation (19) for shock acceleration. Here 

Q(p>9>t)/v Is the number of particles released per (energy-sec) at 
momentum p, longitude 9, time t, and radial distance r-Q from the Sun, and V Is 
the solar wind velocity. 

It Is clear that Fq can be most reliably recovered from f(p,r,t) when Fq 
and the Interplanetary parameters k and V are as uncomplicated as possible. 
Therefore, we consider only events In which the effects of delayed and 
longitude-dependent release at the Sun, longitudinal transport In 
Interplanetary space and all Interplanetary propagation complexities (e.g. 
shocks) are minimal. These are the magnetically well-connected and prompt 
events at solar longitudes In the range 60° ± 40° west of central meridian. 

In addition, the spectrum can be derived with confidence only for those events 
In which the particle fluxes have weak anisotropies during and after the time 
of maximum (Indicating sudden release) and which have smooth time profiles 
(Indicating that the Injection Is uniform over the connection longitudes). In 
such simple events, the particle propagation can be treated as spherical 
diffusion along each field line. 

When, furthermore, k<<Vr and for times t<< r/V, the convection and 
adiabatic deceleration can be neglected. Then equation (30) reduces to simple 
time-dependent spherical diffusion from a point source. If the outer boundary 
Is far away, the solution of equation (30) for k Independent of r Is (Parker 
1965) 
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f (p,r,t) 


Nq(p) exp(-r^/4tk(p)) 
4np^ 2A (tk(p)]'*''^ 


(31) 


where NQ(p)dp is the number of particles between p and p+dp released at the 
Sun. Clearly, If the assumptions so far are valid, f(t) will have the time 
dependence In equation (31), and a fit of particle data to that form will 

confirm the asumptlons and determine Nq(p)* 

A more direct approach Is to note that equation (31) has a maximum at t,^ 

* ^^/6k(p), and that for particles at their tm(^'»p) 
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This Is the basis for the assertion that the TOM spectrum Is the same as that 
released at the Sun. If k « r*^ with b<2, a similar relation holds and N^ip) 
4itp2f(p,r,t„,)r3. 

Many of the assumptions leading to the validity of the TOM spectrum break 
down at low Ion energies because of the long time to maximum (>• 1 day at 1 
MeV/nucleon) and because of the stronger coupling to the solar wind, which Is 
highly variable on time scales of 1 day. In addition, the condition Vr<<k, 
which allows convection and adiabatic deceleration to be neglected, Is 
equivalent to v>>3Vr/X. With X 0.1 AU as In Figure 11 (Palmer 1982), this 
velocity Is '■ 10^ cm/sec, or energies » 1 MeV/nucleon. 

The effect of adiabatic deceleration on the TOM spectrum can be estimated 
several ways, which give similar results. Scholer (1976) has shown by 



54 


numerical simulations of particle trajectories that when k < Vr and k depends 
on energy, convection and adiabatic deceleration not only make the time of 
maximum flux of particles observed at energy E earlier than r2/6k(E), but make 
the discrepancy larger at lower energies. This Is because particles observed 
at energy E had higher energy near the Sun and so diffused faster there; the 
effect Increases at lower energies because the relative energy change Is 
larger at lower energies. The .nergy change must be less than the energy 
shift In the compared to r2/6k(E). This behavior of t^, 1s 

discussed by Kurt et al . (1981) who show that while the tj^ vs. energy curves 
for well-connected flares vary « l/v (as expected for x - constant) above an 
energy Ei (' 1 MeV), the variation Is weaker below Ei In the manner of 
Scholer's (1976) simulations. The break energy E^ varies slightly from flare 

to flare, but Is always near where Vr » k, on the basis to the t^ *^l9h 
energies. Kurt et al . (1981) also show that the TOM spectra of these flares 

are flatter below Ei than above Ei, and they discuss the apparent energy 
change, and Its effect on the TOM spectrum. 

Scholer (1976) found the energy change for each particle by Integrating 
the adiabatic deceleration law dp/dt = - 2Vp/3r along the computed spatial 
trajectories of the particles. Kurt et al. (1981) approximate the 
trajectories with diffusion speed U(p,r) such that dr * Udt. Tnen with U * 
3ak/r, where k « E*, the deceleration law can be Integrated. It can be shown 
that 'J » 2 k/r Is a better choice, consistent with certain analytic results on 
the effect of convection and adiabatic deceleration on the time of maximum and 
the particle flux at that time (Fisk and Axford 1968, Forman 1970j The 
energy change law then becomes dp/dr * - Vp/3k. When x Is constant, this 
results In a change In the speed of a non-relativistic particle, given by Av = 
- vr/x between the Sun and r. For the typical value of X » 0.1 AU, this Is 
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Av " (- 4 x 108 r(AU)) cm/sec. 

Figure 12 shows the effect of this adiabatic energy change on the 
spectrum. 

We conclude that the TOM spectra of solar flare particle fluxes in 
interplanetary space is a good first-order estimate of the spectrum released 
at the Sun, for reasonably well-connected flares showing fairly rapid particle 
release, down to energies where convection and adiabatic deceleration are 
significant. The TOM spectrum should be used because it eliminates the 
velocity dispersion caused by the energy-dependence of the diffusion 
coefficient; the spectrum varies slowly at that time; it is relatively 
insensitive to finite injection times; and it is early enough that convection 
and adiabatic deceleration effects are minimal. The TOM spectrum begins to 
deviate from the spectrum injected at the Sun at energies where k(E) < Vr; 
that is E < El ' 3(r/10\)2 MeV/nucleon for ions and < 1.5 (r/lOx) ^ keV for 
electrons, or when the time of maximum is longer than the solar wind transit 
time divided by the Compton-Getting factor. Typical values of X make ^i ^ 1 
MeV/nucleon for ions and '• 1 keV for electrons. 

The effect of adiabatic deceleration and convection on ion spectra with 
realistic propagation parameters below a few MeV/nucleon is not known 
exactly. We have presented some rules which describe approximately how the 
TOM spectrum is related to the spectrum released at the Sun. Since, as we 
have shown in Section III, the spectrum below 1 MeV/nucleon differs greatly 
between acceleration theories which explain the spectrum above 1 MeV/nucleon 
equally well, a more exact method of deducing the accelerated spectrum below 1 
MeV/nucleon from that observed in interplanetary space is needed to decide 
among such theories. However, because the propagation of such low-energy ions 
is extremely sensitive to interplanetary conditions, it is problematic whether 
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this can be practically done with observations at 1 All. 


V. SUMMARY AND OUTLOOK 


We have reviewed observations and theories of particle acceleration in 
solar flares. Flare accelerated particles that remain trapped at the Sun 
produce a variety of electromagnetic emissions which provide important 
information on the timing of the acceleration and on the number and energy of 
the accelerated particles. While it is clear that not all particles species 
are accelerated to their final energies at the same time, the distinction 
between first-phase and second -phase acceleration is no longer as clear as it 
had been prior to recent gamma-ray observations. There is ample evidence for 
the bulk energization of non-relativistic electrons whose interactions in the 
solar atmosphere produce hard X rays. The energy contained in these electrons 
constitutes a major fraction of all the available flare energy. The 
acceleration of the gamma-ray producing ions in some flares is closely 
associated in time with this bulk energization. The energy content of these 
ions, while smaller than that of the non-relativistic electrons, also 
constitutes a significant fraction of the total flare energy (Section Ila). 

The ions observed in interplanetary space are probably accelerated by 
shocks or turbulence. The correlation between Type II radio bursts and ions 
events provides the main observational support for this conclusion. Both 
stochastic and diffusive shock acceleration can explain the observed energy 
spectra, but the fact that these spectra do not vary much from one event to 
another (Section Ilb) seems to imply that one of these mechanisms is 
dominant. The relative importance of stochastic and shock acceleration 
depends on the injection thresholds, the relative magnitudes of the shock and 
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turbulent speeds and on the magnitudes of the scattering mean free paths. The 
existence of suitable turbulence ahead of the shock is an essential 
requirement for shock acceleration (Section II Ib). 

The chemical compositions of the energetic particles in space and their 
variability appears to be a problem of injection (Section lib). This follows 
from the fact that abundance variations are generally not accompanied by 
spectral variations. Both stochastic and shock acceleration mechanisms have 
injection thresholds, but the threshold of shock acceleration is expected to 
be higher than that of stochastic acceleration (Section III). The abundance 
variations could be due to variations of the threshold energies, or of the 
preheating conditions. Preheating or preacceleration are required to produce 
particles above the thresholds. 

Both stochastic and shock acceleration can be rapid enough to produce the 
observed rise times of the gamma-ray emission. We find that the necessary 
short scattering mean free paths are not inconsistent with other data. 

The question of whether the ions responsible for gamma-ray production and 
those observed in interplanetary space are produced by the same mechanism 
remains unanswered. On the one hand, the energy spectra of these two 
populations appear to be similar, but on the other, the correlation between 
the number of particles seen in space and inferred at the Sun is quite poor. 

An important recent observation (Section Ila) has been that of relativistic 
electrons seen in interplanetary space in correlation with gamma-ray 
emissions. Since stochastic and shock acceleration are not expected to 
accelerate relativistic electron efficiently, this might be <in indication for 
the acceleration of nuclei and electrons to essentially the same energies in 
electric fields (Section IIIc). 

Both observational and theoretical v^ork is required for future progress 



58 


In the understanding of solar flare particle acceleration. On the 
observational side, much progress has recently been achieved by such novel 
Investigations as gamma-ray and neutron observations, X-ra^y observations with 
high temporal resolution, and direct charged particle detections with good 
spectral elemental. Isotopic and Ionic resolutions. These should be continued 
and correlated with each other as well as with other data on solar flares such 
as shock parameters. For example. In energetic particle measurements there Is 
a need for greater sensitivity. Some of the most bizarre and Interesting 
flare events are the small ones; however, existing Instrumentation cannot 
accumulate sufficient statistics to measure the composition In great detail. 
More sensitive gamma-ray observations are needed 1n order to determine whether 
prompt Ion acceleration Is present in all flares. Gamma-ray observations with 
high energy resolution are also needed to resolve the lines and measure 
Doppler shifts which would provide Information on energetic particle 
beaming. 

The theoretical and interpretative work making use of the existing 
observations Is just beginning. We can foresee Important advances from 
studies that will apply recent theoretical results on shock acceleration to 
solar flare acceleration. Investigations of ion acceleration by electric 
fields leading to the beaming of the energetic particles which could have 
observable effects on gamma-ray spectra, and from a much more detailed study 
of acceleration theories which would consider the acceleration of particles 
and their effect on the accelerating agent In a self consistent manner. 
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TABLE 1 


Flare 

Wh(> 1 MeV/n) 

West:(>l MeV/n) 

We(>20 keV) 

No(>10 MeV) 

Nesc.pOlOMeV) 

Aug 4, 1980 

2.5xl030erg 

3xl031erg 

2xl032erg 

1.3x1034 

3x10^5 

June 7, 1972 

2xl029erg 

2xl027erg 

4xl030erg 

8.5x1032 

1031 
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FIGURE CAPTIONS 


Fig. 1. Time dependences of X-rays, gamma-rays and microwave radiation from 
the 4 August 1972 flare (from Bai and Ramaty 1976). 

Fig. 2. Bessel function (see Section Ilia) in velocity, and exponential in 
rigidity fits to proton and alpha particle spectra for the 5 November 1974 
flare particle event. Fit parameters are aT (proton) = 0.024, aT (alpha) = 
0.015, Rq (proton) = 73 MV, Rq (alpha) = 80 MV. 

Fig. 3. Electron energy spectrum from the 7 September 1973 flare (from Lin, 
Mewaldt and Van Hollebeke 1982). 

Fig. 4. Correlation between solar flare electron and proton intensities 
observed in interplanetary space (from Ramaty et al. 1980). 

Fig. 5. Compositon of solar energetic particles relative to the photosphere 
(see Section Ilb for details). 

Fig. 6. Energy loss rates of protons, Fe nuclei and electrons in neutral H 
(dashed curves) and energy gain rates from equation (18) (from Ramaty 1979). 

Fig. 7. Energy loss rates of protons and electrons in an H plasma of 2x1 0&K 
(dashed curves). The loss rates for other ions scale as Z2 /a. The energy 
gain rates are from equation (18) (from Ramaty 1979). 

Fig. 8. Time dependent proton and electron spectra from stochastic 
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acceleration with impulsive injection and no losses, Eq injection 

energy (from Ramaty 1979). 

Fig. 9. Comparison of the particle spectra produced by steady-state stochastic 
acceleration, equation (13), with shock acceleration, equation (24), for 
models in which the scattering mean free path, x, and the escape time, T, 
are energy independent. Vgh the shock speed and (5V)2 the mean square 
turbulent velocity. 

Fig. 10. Magnetic field and flow configuration in (a) the Petchek (1964) and 
(b) the Sonnerup (1970) reconnection models. 

Fig. 11. Mean free paths of solar electrons and protons in interplanetary 
space deduced by fitting time profiles in many events (from Palmer 1982). 

Fig. 12. Illustration of the probable effect of adiabatic deceleration on the 
spectra of solar flare ions at the time of maximum. The spectrum at the Sun 
(shown as a dashed line) is shifted to lower energies. 
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